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Abstract. Mixing Mg with Ti leads to a hydride Mg^Ti(i_3,)Il2 with markedly 
improved (de)hydrogenation properties for x < 0.8, as compared to MgH2. Optically, 
thin films of Mg^Ti(i_2,)H2 have a black appearance, which is remarkable for a hydride 
material. In this paper we study the structure and stability of Mg^Ti(i_2,)H2, x = 0- 
1 by first-principles calculations at the level of density functional theory. We give 
evidence for a fluorite to rutile phase transition at a critical composition Xc — 0.8-0.9, 
which correlates with the experimentally observed sharp decrease in (de)hydrogenation 
rates at this composition. The densities of states of Mg^Ti(i„3,)Il2 have a peak at 
the Fermi level, composed of Ti d states. Disorder in the positions of the Ti atoms 
easily destroys the metallic plasma, however, which suppresses the optical reflection. 
Interband transitions result in a featureless optical absorption over a large energy 
range, causing the black appearance of Mg^Ti(i_2,)Il2. 
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1. Introduction 

Hydrogen promises to be a good candidate to replace carbon based fuels in the future [1] . 
The current inability to store hydrogen in a safe and sufficiently dense form obstructs 
its use in vehicles or in portable applications. A practical means of storage should yield 
a system that operates at moderate temperatures and releases hydrogen fast enough to 
feed a fuel cell [2] . Storage should be reversible and the dehydrogenated system should 
adsorb hydrogen sufficiently fast. Typical target numbers for gravimetric and volumetric 
hydrogen densities in storage systems are 6 weight % (wt%) and 0.045 kg/L. 

The highest volumetric hydrogen densities among all possible storage forms are 
achieved in metal hydrides, where atomic hydrogen is bonded in the crystal lattice of 
bulk metals or alloys [31 H]. Reasonable gravimetric hydrogen densities can be obtained if 
lightweight metals are used. For instance, the simple dihydride MgH2 has a gravimetric 
hydrogen density of 7.7 wt%. However, MgH2 is so stable that releasing hydrogen at a 
pressure of 1 bar requires a temperature of around 300°C [5], which is impractically high 
for use in combination with PEM fuel cells. Moreover, the (de)hydrogenation kinetics 
of MgH2 is poor. Surface processes can hamper the kinetics, as the dissociation rate 
of hydrogen molecules on a Mg surface is small [6] , and surface oxidation leads to the 
formation of a hydrogen diffusion barrier [7]. The surface kinetics can be improved by 
annealing [8J and by adding a catalyst 0. Experimental studies also claim that the 
(de)hydrogenation kinetics of bulk MgH2 is very slow. It has been suggested that the 
kinetics is hampered by the rutile crystal structure of a;-MgH2 [IDl [HI [12] . 

It is well known that 3d transition metals (TMs) act as catalysts to improve the 
surface kinetics of hydrogen adsorption in Mg [13]. To obtain this improvement usually 
only a few wt% of TM is added to Mg. Recently it has been shown that also the bulk 
(de)hydrogenation kinetics can be improved markedly by adding TMs. One has to add 
a substantially larger amount, however, and make alloys Mg^.TM(i_a;) with TM = Sc, 
Ti and x < 0.8 [T0[ [Til [I2]- The hydrides of these alloys have a cubic crystal structure 
quite unlike the a-MgH2 rutile structure [m [151 [13 [IZ]. Obviously, to preserve a 
high gravimetric hydrogen density it is essential to use lightweight TMs. As Sc is too 
expensive to be used on a large scale, Ti is then the obvious choice. Thin Mg2.Ti(i_a,) 
films, prepared by various experimental techniques, can be reversibly hydrogenated and 
dehydrogenated [TH UHl |19l [20]. The (de)hydrogenation kinetics is fast and increases 
gradually for alloys with increasing Mg content x, until a maximum at x ~ 0.8. If the 
Mg content is larger than this value, the kinetics becomes much slower. 

Recent experimental studies reveal that Mg-Ti-H thin films have interesting optical 
properties, which may lead to applications of these compounds other than hydrogen 
storage [T9|, [201 [21] . In the dehydrogenated state the films are highly refiective, whereas 
upon hydrogenation they become black, i.e. they have a low refiection and high 
absorption for light in the visible range. This would allow them to be used as switchable 
smart coatings for solar panels, or as hydrogen sensors to detect the presence of hydrogen 
gas, for instance. 
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Despite the experimental interest in Mg-Ti hydrides, obtaining structural data 
has proven to be difficult as so far the bulk compounds have resisted synthesis under 
normal conditions. The only bulk Mg-Ti-H compound synthesized so far is MgyTiHig, 
obtained by letting the dihydrides MgH2 and TiH2 react under extreme conditions 
of high pressure (8 GPa) and temperature (873 K) [221 [23] . In contrast, thin films 
of Mg^.Ti(i_a,)H2 can be readily obtained for any composition x by various deposition 
techniques. The limited amount of data on the structure, the thermodynamic stability 
and the electronic properties of Mg-Ti-H compounds has motivated the present first- 
principles computational study on Mg^Ti(i_a;)H2. 

First we consider the simple dihydrides MgH2 and TiH2 and subsequently 
characterize the properties of Mg^Ti(i_^.)H2, focusing on the relative stability of possible 
cubic and rutile crystal structures. Limited first-principles calculations on ordered 
structures have indicated that a cubic fiuorite-like crystal structure of Mg^TM(i_2.)H2 is 
more stable than the rutile structure for a range of early TMs and x < 0.8 [2H[25]. Here 
we focus on TM = Ti and study such structures in more detail, including the effects 
of disorder in the positions of the metal atoms, which is relevant for the experimental 
thin film studies. Our results suggest that the observed change in kinetics as a function 
of Mg content, is associated with a change in crystal structure. Finally, we calculate 
the optical properties of Mg3.Ti(i_j.)H2 and analyze them in terms of the compounds' 
structure and electronic structure. We show that the optical black state, which makes 
these compounds unique among the reversible hydrides so far, is an intrinsic property 
[26] . unlike similar optical states of other Mg-TM-H thin films, which are caused by 
multiphase inhomogeneities in thin films [27t [28]. 

2. Computational Details 

Most of the present calculations are performed at the level of the generalized gradient 
approximation (GGA) to density functional theory, using the projector augmented wave 
(PAW) technique and a plane wave basis set [291 [30], cis implemented in the Vienna Ab 
initio Simulation Package (VASP) [31], [32|, [33] . Standard frozen core PAW potentials 
are used and the H Is, Mg 2s, Ti 4s and 3d shells are treated as valence shells. As Ti 
has a partially filled 3d shell we have considered the possibility of spin polarization and 
ferromagnetic or antiferromagnetic orderings. Unless otherwise mentioned explicitly, 
most of the compounds studied turn out to be paramagnetic, however. The main 
results are obtained using the Perdew Wang 91 (PW91) functional [3l], but we have also 
performed some tests with the hybrid DFT/Hartree-Fock schemes B3LYP and HSE06 

[351 [Ml [ag [an [39]. 

We require that the total energy of each compound is converged to within 1 
meV/atom with respect to the plane wave kinetic energy cutoff, which is assured by 
using a cutoff of 650 eV. Convergence with respect to the fc-point sampling for the 
Brillouin zone (BZ) integration is tested independently on simple compounds using 
regular meshes of increasing density. We aim at converging total energies on a scale of 1 
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meV/atom. This is obtained with the following /c-point meshes for the simple unit cells: 
12x12x8 for hep unit cells, 16x16x16 for cubic, fee or fluorite, unit cells and 12x12x16 
for bet or rutile unit cells. In studying alloys and their hydrides we construct supercells 
of these simple unit cells and keep the same fc-point grid. Thus employing the equivalent 
k-points method one avoids relative /^-point sampling errors and maintains the accuracy 
for all cell sizes |40j . 

The molar volume of Ti (10.64 cm^) is much smaller than that of Mg (14.00 cm^), 
which makes finding the optimal cell parameters for Mg^Ti(i_a;)H2 non-trivial. We start 
by keeping the atomic fractional coordinates and the c/a ratio fixed and first carry out 
volume scans for each composition x and crystal type individually. Values for the lattice 
parameters are then estimated from the total energy minima. These are then used as 
starting values for a full optimization of each composition and crystal type. 

The cell parameters, including the cell volume, and the atomic positions within the 
cells are optimized by minimizing the forces and stresses with the conjugate gradient 
algorithm [H] . During optimization Methfessel-Paxton smearing is used with a smearing 
parameter of 0.1 eV [42]. The criteria for self-consistency are set to 10^^ eV and 10~^ 
eV for the total energy differences between two consecutive electronic and ionic steps, 
respectively. Structural relaxation is assumed to be complete if all the forces acting 
on atoms are smaller than 1 meV/A in the simple cell studies, and 10 meV/A in the 
supercell studies. After the structures are optimized, the total energies are recalculated 
self-consistently with the tetrahedron method [43]. The latter technique is also used to 
calculate the electronic density of states (DOS). 

The properties of the H2 molecule are calculated in a cubic cell with cell parameter 
10 A, using F-point sampling. The calculated H— H bond length, binding energy, and 
vibrational frequency are 0.748 A, —4.56 eV and 4351 cm~^, respectively, in good 
agreement with the experimental values of 0.741 A, —4.48 eV and 4401 cm~^ [69]. 

To study the stability of Mg^Ti(i_^)H2 we look at the reaction 

xMg + (1 - x)Ti + H2(gas) ^ Mg,Ti(i_,,)H2 (1) 

In principle one should consider the change in Gibbs free energy G{T, P) = U + PV — TS 
of this reaction. As the thermodynamic properties of hydrogen gas are well documented 
|44j . one can focus upon the solid phases involved in the reaction. In addition, 
the temperature dependence of G of such solids tends to be relatively small over 
the temperature range of interest [lH HH]. A calculation of the change in enthalpy 
H = U + PV at T = OK should therefore be adequate to assess the relative stability 
of metal hydrides. Since moreover the PV contribution of solids can be neglected, it 
suffices to focus upon the change in energy U. 

U corresponds to the total energy of a system as obtained in a DFT calculation. 
The positions of the atomic nuclei are then fixed, which corresponds to T = K for 
classical particles. However, hydrogen is such a light atom that quantum vibrational 
energies give contributions to the total energies that are not negligible, even at T = 
K. These zero-point energies (ZPEs) associated with atomic vibrations are obtained 
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by solving the eigenvalue problem of the dynamical matrix [17]. Slightly displacing an 
atom from its equilibrium position creates forces on the other atoms. The dynamical 
matrix can then be constructed by displacing each atom at a time and finite differencing 
the forces. We use symmetric displacements dr = ±0.01 A for each atom and its three 
degrees of freedom. These calculations are carried out in supercells to properly account 
for the spatial range of the dynamical matrix. 

The ZPE corrected reaction energies are then calculated as 

p r 

where r,p indicate the reactants and products of the reaction, respectively, U is the 
total energy obtained from a DFT calculation and f/^'*^ is the vibrational ZPE. As it 
turns out, the ZPE corrections result in a rather constant shift of the reaction energy 
for the range of compounds studied here. Therefore, we frequently give the reaction 
energy as AU, calculated without Up^^ contributions, and mention the ZPE corrections 
separately. 

As we will discuss in Sec. 14. we employ two different strategies to model the 
structures. The first strategy consists of constructing a set of relatively simple ordered 
structures, inspired by the experimental structure of MgyTiHie [22l [23] . The structure of 
Mg3;Ti(i_^)Hj, as deposited in thin films is cubic with Mg and Ti atoms at fee positions, 
but without a regular ordering of Mg and Ti atoms at these positions [T5[[T7]. To model 
such disordered structures, our second strategy consists of employing special quasi- 
random structures (SQSs), which enable to model random alloys in a finite supercell 
[I8] . At each composition Mg^Ti(i_a;) of interest, the Mg and Ti atoms are distributed 
such, that their lower order correlation functions are equal to those of a perfect random 
alloy [50] . Hydrogen atoms are then inserted into similar positions as in the simple 
structures. We use a set of SQSs for each crystal structure (fluorite or rutile) and 
composition, with a total number of atoms ranging from 48 to 192. Again all cell 
parameters are optimized, as well as the positions of all atoms within the cell. 

To model the optical response of Mg2.Ti(i_^)H2 we calculate the dielectric function, 
consisting of interband and intraband contributions 

e{uj) = Einteriuj) + £:mtra(t^), (3) 

which are calculated separately. The interband contribution eintor is calculated in the 
independent particle approximation, using DFT eigenvalues and wave functions, i.e. 
neglecting exciton and local field effects. Im[einter(i^)] is calculated directly via the 
standard longitudinal expression, and Re[£:intcr(^)] is obtained by a Kramers-Kronig 
transform [5T1 [52l [26] . As optical data on hydrides are usually obtained from micro- 
and nano-crystalline samples whose crystallites have a significant spread in orientation, 
we use the directionally averaged dielectric function. 

Many metals and metal alloys undergo a metal-insulator transition upon full 
hydrogenation [53l[27l[28l[Mll55ll56l[52l[57l. Mg^Ti(i_^)H2 however retains a finite DOS 
at the Fermi level. Hence intraband transitions contribute to the dielectric function. We 




Figure 1. Crystal structures of Q;-MgH2 (a), /3-MgH2 (b), and /-TiH2 (c). Mg, Ti 
and H atoms are shown as green, blue and white coloured spheres, respectively; the 
polyhedra indicate the coordination of the metal atoms by hydrogens. The images are 
made using the VESTA software [HT| . 



model Eintrai^) by a standard free electron plasma model, with the plasma frequency 
calculated from the k-derivatives of the DFT eigenvalues [58ll26]. Once the full dielectric 
function e{uj) is calculated, optical functions such as the extinction coefficient or 
the refractive index n{uj) are obtained via the standard expression 

e{uj) = [n{uj) +iK{uj)f (4) 



3. Simple dihydrides 

As discussed in the introduction, experimental information on the structure of 
Mg2.Ti(i_j.)H2 is limited, in particular concerning the positions of the hydrogen atoms. 
A sensible way to approach the problem is to assume that the structure is similar either 
to MgH2 or to TiH2. In this section we consider these simple dihydrides in more detail. 

At standard temperature and pressure MgH2 has the rutile structure, space group 
P42/mnm (136), see figure [H In this so-called a-MgH2 phase, six hydrogen atoms 
surround each Mg atom in a distorted octahedron. The distance between Mg and two 
out of the six surrounding H atoms is 1.938 A, and between Mg and the other four H 
atoms it is 1.951 A. Each H atom is in the centre of a triangle with Mg atoms at the 
vertices. Experimentally it has been shown that an increase in the temperature and 
pressure to 923 K and 4 GPa, converts a-MgH2 to a cubic form, called /5-MgH2, with 
no further information about the positions of the H atoms [59]. DFT calculations have 
suggested a structure for /3-MgH2 [SO] • Starting from a fiuorite cubic structure with the 
hydrogen atoms on tetrahedral interstitial sites, space group Fm3m (225), structural 
optimization gives a cubic structure with reduced symmetry, space group Pa3 (205). 
This is accompanied by a decrease in the total energy of 0.24 eV/f.u. and a volume 
expansion of ~ 4%. Each hydrogen atom is displaced from the centre of a tetrahedron 
to a base plane, where the coordination by Mg atoms is triangular with the three Mg 
atoms at a distance di = 1.911 A, similar to the a-phase. The fourth Mg atom of 
the original tetrahedron is at a distance d2 = 2.878 A. Each Mg atom is octahedrally 
coordinated by six H atoms at a distance di. An additional two H atoms are at d2- At 
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Table 1. Calculated and experimental structural parameters of simple metals and 
hydrides. 

Space group 



(/DiTi "nm inn 


Unit cell (A) 








V 

J 


2 


Mg 


PGa/mmc (194) 


Mg 


2c 


1/3 


2/3 


1/4 




a ^ 3.190 ; c 5.180 












exp.[7Dj 


a = 3.209 ; c = 5.210 


Mg 


2c 


1/3 


2/3 


1/4 


Ti 


PGa/mmc (194) 


Ti 


2c 


1/3 


2/3 


1/4 




a = 2.916 ; c = 4.631 












exp.[70] 


a = 2.951 ; c = 4.686 


Ti 


2c 


1/3 


2/3 


1/4 


a-MgH2 


P42/mnm (136) 


Mg 


2a 













a — 4.4y4 , c — o.uuo 


TT 

rl 


/If 


U.OU44 


U.oU44 


U 


exp.[6I] 


a = 4.501 ; c = 3.010 


Mg 


2a 















H 


4f 


0.3044 


0.3044 





/3-MgH2 


Pa3 (205) 


Mg 


4a 













a = 4.796 


H 


8c 


0.3464 


0.3464 


0.3464 


/-TiH2 


Fm3m (225) 


Ti 


4a 













a = 4.424 


H 


8c 


1/4 


1/4 


1/4 


exp.[70] 


a = 4.454 


Ti 


4a 















H 


8c 


1/4 


1/4 


1/4 


r-TiH2 


P42/mnm (136) 


Ti 


2a 













a = 3.991 ; c = 3.080 


H 


4f 


0.322 


0.322 






zero pressure the optimized /3-MgH2 structure is energetically unfavourable compared to 
the Q;-MgH2 structure by 0.10 eV/f.u. and its volume is ~ 9% smaller. The optimized 
structural parameters obtained in our calculations are given in table [H They are in 
good agreement with available experimental data. 

The most stable phase of TiH2 at standard temperature and pressure has the fluorite 
structure, space group Fm3m (225), see figured! where the Ti atoms occupy the fee 
lattice positions and hydrogen atoms reside at the tetrahedral interstitial positions. 
The calculated value of the Ti— H bond length is 1.916 A. We call this the /-TiH2 
structure and to have a comparison between the two structural types of dihydrides we 
also construct a rutile form of TiH2, and call it r-TiH2 in table [H This artificial rutile 
structure does not exist in nature to our knowledge. It has a 13% larger volume than 
/-TiH2, and its total energy is 0.65 eV/f.u. higher. 

The formation energies of the simple hydrides are defined as the reaction energies 
corresponding to ([T]) for x = 0, 1. The calculated values with and without vibrational 
ZPE corrections are given in tableO ZPEs tend to destabilize the metal hydrides. By far 
the biggest contribution stems from the difference between the vibrational frequencies 
of the hydrogen atoms in the metal hydride and in the gas phase H2 molecules. The 
ZPEs of the metals are small, i.e. 0.031 and 0.033 eV/atom for Mg and Ti, respectively 
(calculated with a 3x3x2 hep supercell). In general, the more compact the metal 
hydride structure, the less space for the hydrogen atoms to vibrate, and the higher 
their vibrational frequencies and the ZPE are. The volumes per formula unit are 30.34 
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Table 2. Calculated formation energies in eV without (A?7) and with (A[/^^^) 
vibrational zero point energy corrections, cf. ([2]), compared to experimental values 
(A[/°''P). The latter are taken from [71] unless indicated otherwise. 



Compound 


AU 








a-MgH2 


-0.657 


-0.558 


-0.682 


-0.694 -0.771 








-0.772 


[M] -0.779 [zg 


a-TiH2 


-1.467 


-1.279 


-1.297 


-1.351 -1.378 








-1.450 


[73^ -1.496 



Table 3. Calculated formation energies of Q;-MgH2 in eV, obtained with different PAW 
pseudopotentials for Mg and different DFT functionals or hybrid DFT/Hartree-Fock 
schemes. The ZPE correction has been calculated at the PW91 level only. 



Mg valence 


functional 


AU 


AC/2PE 


3s 


PW91 


-0.657 


-0.558 


2p3s 


PW91 


-0.648 


-0.549 


3s 


HSE06 


-0.25 


-0.15 


3s 


B3LYP 


-0.95 


-0.85 



for a-MgHa and 21.64 for /-TiHg. The calculated ZPEs are 0.406 eV/f.u. and 
0.501 eV/f.u., respectively. These are calculated using a 2x2x3 rutile supercell for 
MgH2, containing 72 atoms, and a 2x2x2 fluorite supercell for TiH2, containing 96 
atoms. The calculated ZPE of a H2 molecule is 0.276 eV. Consequently, including ZPE 
corrections destabilizes a-MgH2 by approximately 0.1 eV, whereas the more compact 
/-TiH2 is destabilized by almost 0.2 eV, see table [2l 

The experimental numbers obtained for the formation energy of TiH2 range from 
— 1.30 to —1.50 eV. The calculated numbers are in good agreement with these results. 
Experimental values of the formation energy of MgH2 range from —0.68 to —0.78 eV. 
The calculated values overestimate these numbers somewhat, in particular if one includes 
ZPEs. Such differences between calculated and experimental formation energies are also 
found in other simple alkali and alkaline earth hydrides and seem to be typical for the 
use of GGA functionals such as PW91 or PBE. To see whether the calculated values 
can be improved, we test the use of a hard PAW pseudopotential for Mg, which includes 
the 2p shell as valence electrons, as well as the use of the hybrid DFT/Hartree-Fock 
schemes B3LYP and HSE06. The results, listed in table [3|, show that the effect of using 
a hard PAW potential is minor. Hybrid DFT/Hartree-Fock schemes do not necessarily 
give an improvement over DFT/GGA functionals. The HSE06 and B3LYP schemes lead 
to formation energies of MgH2 that are respectively much larger and somewhat smaller 
than the experimental one. Since in the following we are foremost interested in energy 
differences between crystal structures, we stick to using the PW91 GGA functional. 

a-MgH2 is an insulator with an experimental optical band gap of 5.6 eV [62]. As 
usual DFT calculations underestimate this gap [52|; the value obtained for Q;-MgH2 is 
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Figure 2. Electronic density of states (DOS) of /3-MgH2 and /-TiH2. The origin of 
the energy scale is set at the top of the valence band. The top figure gives the total 
DOS, and the middle and bottom figures represent the local DOS, projected on the 
atoms; s, p and d contributions are shown in red, green and blue, respectively. The 
regions under the s curves are shaded and the d curves are dashed. 

3.8 eV, whereas the calculated gap for /?-MgH2 is 2.7 eV. The electronic (projected) 
densities of states ((P)DOS) of the cubic simple dihydrides /5-MgH2 and /-TiH2 are 
shown in figure El The fact that many fully hydrided simple metals and alloys are 
insulators can be understood from simple chemical arguments [SH [551 EHl [521 [SI!- The 
Mg atoms almost fully donate their valence (2s) electrons to the H atoms. Indeed the 
PDOS in figure [2] shows that the valence bands have a predominant H Is character, 
whereas the Mg 2s states emerge in the conduction bands. In contrast, TiH2 is metallic 
with a considerable DOS at the Fermi level. The PDOS shows that the states around 
the Fermi energy have a dominant Ti 3d character, whereas, as for MgH2, the lower 
valence bands are dominantly H Is. It suggests that in TiH2 two electrons per metal 
atom are transferred to hydrogen and, as Ti atoms have four valence electrons (s^d^), 
they remain in an open-shell configuration in the hydride. Figure [2] illustrates that the 
Ti d electrons participate significantly less than the s electrons to the bonding to H 
atoms. 

4. Mg,Ti(i_,)H2 
4.I. Structures 

The only detailed structure of Mg^Ti(i_^)H2 available from experiment is for the 
composition x = 0.875. In the MgyTiHig high pressure phase the metal atoms are 
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Table 4. Optimized cell parameters and atomic positions of the fluorite-type simple 
Mg^Ti(i_j,-)H2 structures. 

Compound, Space group 



Unit cell (A) 






X 


V 

J 


z 


M^So. 125^10. 875H2 


Mg 


4a 











Fm3m (225) 


Til 


4b 


1/2 


1/2 


1/2 


a = 8.8896 


Ti2 


24d 





1/4 


1/4 




HI 


32/ 


0.1285 


0.1285 


0.1285 




H2 


32/ 


0.3755 


0.3755 


0.3755 


Mgo.25Tio.75H2 


Ti 


3c 





1/2 


1/2 


Pm3m (221) 


Mg 


la 











a = 4.4648 


H 


85 


0.2566 


0.2566 


0.2566 


Mgo.5Tio.5H2 


Til 


la 











P4/mmm (123) 


Ti2 


Ic 


1/2 


1/2 





a = 4.4720 


Mg 


2e 





1/2 


1/2 


c = 4.6544 


H 


8r 


0.25 


0.25 


0.2352 


Mgo.75Tio.25H2 


Ti 


la 











Pm3m (221) 


Mg 


3c 





1/2 


1/2 


a = 4.6208 


H 


85 


0.2432 


0.2432 


0.2432 


Mgo.875Tio.l25H2 


Ti 


4a 











Fm3m (225) 


Mgl 


4& 


1/2 


1/2 


1/2 


a = 9.3565 


Mg2 


24d 





1/4 


1/4 




HI 


32/ 


0.1194 


0.1194 


0.1194 




H2 


32/ 


0.3725 


0.3725 


0.3725 


Mgo.875Tio.i25H2 exp.[23] 


HI 


32/ 


0.094(2) 


0.094(2) 


0.094(2) 


a = 9.564(2) 


H2 


32/ 


0.365(2) 


0.365(2) 


0.365(2) 



in fee positions and are ordered as in the CayGe structure [22l [23]. The H atoms are at 
interstitial sites, but displaced from their ideal tetrahedral positions. This structure has 
96 atoms in the unit cell. We use it as a starting point to construct relatively simple, 
fluorite-type structures for other compositions x. For x = 0.125, 0.875 we use the CayGe 
structure to order the metal atoms, for x = 0.25, 0.75 the CusAu (-^^12) structure, and for 
X = 0.5 the CuAu (i^lo) structure. The H atoms are placed at or close to tetrahedral 
interstitial positions. In each of the structures and compositions the cell parameters 
are optimized, as well as the positions of all atoms within the cell. Care is taken to 
allow for breaking the symmetry in the atomic positions, as the hydrogen atoms are 
often displaced from their ideal tetrahedral positions. Although strictly speaking such 
structures are then no longer ideal fiuorite structures anymore, we still use the term 
fluorite in the following. 

The optimized fiuorite structures are given in table HI Most of these structures are 
cubic, except MgQ 5Tio.5H2, which is tetragonally distorted. Comparing the calculated 
and the experimental structures of Mgg g75Tio.i25H2 one observes that the calculated 
lattice parameter is 2% smaller than the experimental one. Although such a result 
is not uncommon for DFT calculations, note that the differences between calculated 
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0.125 0.25 



0.5 
X 



0.75 0.875 1 

P- MgH^ 



Figure 3. The volumes per formula unit (V/f.u) in A'^ of Mg^Ti(i_^)H2 in the ordered 
fluorite structures, as a function of the composition x. The dashed line marks a linear 
interpolation between /-TiH2 and /3-MgH2. The crosses mark the volumes of the 
special quasi-random structures. 



and experimental lattice parameters of MgH2 and TiH2 are much smaller, see table [TJ 
The Wyckoff positions and the hydrogen-metal coordinations of the calculated and 
experimental structures of Mgg 875Tio.i25H2 agree, but the absolute difference between 
the calculated and experimental H positions is 0.1-0.4 A. This is not a volume effect, 
since optimizing the structure using the experimental or the calculated lattice parameter 
gives the same H positions. We have also checked that breaking the symmetry by random 
displacement of the H atoms yields the same optimized Fm3m (225) structure. From 
our calculations the experimental structure is 0.766 eV/f.u. higher in energy. 

The calculated tetrahedral environment of a HI atom in Mgg 875Tio.i25H2 consists 
of a Ti atom at a distance of 1.935 A and three Mg2 atoms at 2.058 A; the experimental 
distances are 1.557 A and 2.293 A, respectively. The calculated tetrahedral environment 
of a H2 atom consists of a Mgl atom at 2.066 A and three Mg2 atoms at 2.013 A 
compared to experimental distances of 2.236 A and 2.021 A. The calculated H— Ti 
distance is close to the 1.916 A found in TiH2, see tabled], whereas the experimental 
distance is significantly smaller. In fact, for all compositions in table H] we find H— Ti 
distances in the range 1.91-1.94 A. From the calculated H— Mg distances one observes 
that the H atoms are closer to ideal tetrahedral positions than in the experimental 
structure. From the positions of the H2 atoms it is clear, however, that both the 
calculated and the experimental structure are closer to an ideal fluorite structure than 
to /3-MgH2. 

One should note that the experimental structure of Mg7TiHi6 is based upon room 
temperature XRD data. As XRD measures the electronic charge distribution, modelling 
the data with atomic scattering functions introduces errors on the positions of light 



First principles modelling of magnesium titanium hydrides 12 




Figure 4. Cells to model special quasi-random structures of Mg^Ti(i_a;)H2 for 
X = 0.75 (left image) and x = 0.875 (right image) compositions. Mg, Ti and H 
atoms are shown as green, blue and white coloured spheres, respectively; the polyhedra 
indicate the coordination of the metal atoms by hydrogens [61) . 

atoms such as hydrogen. In particular, the asymmetry in the hydrogen positions is 
sometimes exaggerated. In that case low temperature neutron diffraction can give a 
structure that is much closer to the calculated structure [631 SS] • 

The calculated volumes of the fiuorite structures of Mg^Ti(i_2,)Il2 as a function of 
the composition x, normalized per formula unit, are shown in figure O According 
to Zen's law of additive volumes one would expect a linear dependence [M]- The 
curves deviate slightly, but distinctly, from straight lines, with a maximum deviation of 
5.4 % at X = 0.5. This deviation is consistent with the experimental observations on 
Mg^Ti(i_2.)H2 [ISlllS]. It is also observed in simple metal alloys [6^ . 

The binary phase diagram of Mg and Ti indicates that they do not form a stable 
bulk alloy. It is however claimed that synthesis by mechanical alloying may be possible 
|65j . In addition, thin films of Mg^Ti(i_a,) are readily made, which can be reversibly 
hydrogenated [IDl [El [HI [ISl [20]. The crystal structure of Mg^Ti(i_^.)Hy in thin 
films, X < 0.8, y ~ 1-2, is cubic, with the Mg and Ti atoms at fee positions, but no 
detectable regular ordering of Mg and Ti atoms at these positions [151 [17]. The positions 
of the H atoms have not been determined from experiment. 

To model such disordered structures, we perform calculations on SQSs, which enable 
to model random alloys in a finite supercell. We use a 32 atom supercell to model SQSs 
of fee Mga.Tii_a; for X = 0.25, 0.5, and 0.75, and a 64 atom supercell for 0.875 [19]. 
Inserting hydrogen atoms in tetrahedral interstitial positions then gives supercells with 
a total number of atoms of 96 and 192, respectively. Optimizing the structure leads to 
some displacements in the positions of the metal atoms, but they remain close to fee. 
There is a larger spread in the positions of the H atoms, which increases with increasing 
X. Representative examples of SQSs are given in figure [H Whereas in the ordered 
structures of table |H the nearest neighbour Ti— Ti distance increases with increasing 
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Figure 5. Radial distribution functions (RDFs) for Mg— H (left column) and Ti— H 
(right column) pairs in Mg^Ti(i„j;)H2, x = 0.75, 0.875 fluorite-type SQSs, compared 
to the RDFs of /3-MgH2 and /-TiHs. The red lines indicate the integrated RDF. 
Horizontal dotted and dashed lines indicate the H coordination numbers of Mg and Ti 
in the simple dihydrides, respectively. 

Mg content x, in the SQSs one can always find at least one Ti— Ti pair at a distance 
comparable to that in /-TiH2, for the compositions studied. 

The optimized volumes of the SQSs are shown in figure [3l For x < 0.75 these 
volumes are within 1.5% of those of the ordered structures. Remarkably, the SQS 
volume of MgQ 875Tio.i25H2 is 4.9% larger than that of the ordered structure of table HI 
This is indicative of a structural change at this composition, which can be clarified 
by calculating the radial distribution functions (RDFs) of the SQSs. The RDFs of 
Mg3,Ti(i_j,)H2 for x = 0.75, 0.875 are plotted in figure [5l Shown are the RDFs of Ti— H 
and Mg— H distances and the corresponding RDFs in the cubic /3-MgH2 and /-TiH2 
compounds. Integrating the RDF gives the number of hydrogen atoms surrounding 
each metal atom. The number in the first shell then corresponds to the coordination 
number of the metal atom. 

In /-Till2 each Ti atom is coordinated by eight H atoms at a distance of 
approximately 1.9 A, whereas in the Mg atoms in /3-MgH2 are coordinated by six H 
atoms at 2.0 A, and by an additional two H atoms at 2.9 A. The RDFs of the SQSs 
show that in Mg^Ti(i_2,)H2 each Ti is coordinated as in /-TiH2, irrespective of the 
composition x. For x = 0.75 each Mg atom is also coordinated by eight atoms at a 
similar distance, as in a fluorite structure. However, for x = 0.875 the RDF for Mg— H 
starts to resembles that of /3-MgH2, i.e. a peak representing six H atoms at ~ 2.0 A, 
and a shoulder representing two H atoms at ~ 2.9 A. Apparently a large Mg content 
is required before such a structure can develop, which is accompanied by a volume 
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expansion with respect to an ideal fluorite structure. In order not to complicate the 
terminology we call all such structures fluorite-type in the following. 

Since we will compare the stability of these fluorite-type structures to rutile-type 
structures Mg^Ti(i_j,)H2, we also make models of the latter using the a-MgH2 structure 
as a starting point. To construct simple, ordered structures, we replace a fraction x of 
the Mg atoms by TM atoms and use the smallest supercell of the rutile structure where 
this leads to an integer number of atoms. In case of multiple possible cells, the results 
given below refer to the cell that leads to the lowest energy. Rutile SQSs of Mg^Ti(i_2,) 
for X = 0.25-0.875 are constructed using standard algorithms [19]. As it turns out that 
rutile structures are relatively unstable over a large composition range, we refrain from 
giving their structural details here. 

4-2. Stability 

Experimentally it is concluded that Mg and Ti do not form a thermodynamically stable 
bulk alloy. Indeed we find that its costs energy to make alloys from the bulk metals 
Mg and Ti. For instance, the energy costs of making Mg2,Ti(i_j;) fee SQSs from bulk 
hep Mg and Ti are 0.16 and 0.11 eV/atom for x = 0.5 and 0.75, respectively, and for 
Mg^Ti(i_a:) in the simple, ordered CuAu and CusAu structures the energy costs are 0.22 
and 0.15 eV/atom, respectively. Calculations on hep Mg^Ti(i_j.) decrease these numbers 
by ~ 0.04 eV, which implies that the hep alloy is more stable than the fee alloy, but it 
is still unstable. Moreover, simple estimates of the configurational entropy show that 
the alloys cannot be stabilized at a reasonable temperature by entropy. Metastable 
structures can still be relevant, though, as thin film deposition of Mg2.Ti(i_i.) yields a 
hep structure that survives a large number of (de)hydrogenation cycles. 

The formation energy of Mg^Ti(i_a;)H2 according to ([1]) is given in figure [6] for the 
fluorite and rutile structures. The numbers are given without ZPE corrections, which we 
will discuss later. Clearly all numbers are negative, which implies that the compounds 
are stable with respect to decomposition into the elements. In addition, over most of 
the composition range the fluorite structure is more stable than the rutile structure. 
For X larger than a critical composition Xc the rutile structure becomes lower in energy. 
If we use the energies calculated for the simple ordered structures as a starting point 
and linearly interpolate between the compositions, we find Xc ~ 0.83. The fact that Xc 
has a high value makes sense, since the energy difference between the Q;-MgH2 (rutile) 
and /3-MgH2 ("fluorite") phases is only 0.10 eV/f.u., whereas between /-TiH2 (fluorite) 
and r-TiH2 (rutile) it is 0.65 eV/f.u.. Therefore, it requires a high Mg content to force 
the structure into rutile. Indeed, experimental results suggest that the cross-over from 
fluorite to rutile structures in Mg^Ti(i_a,)H2 occurs for x somewhere in the range 0.8-0.9 

mm- 

Note that for all compositions the structures are metastable with respect to 
decomposition into Till2 and MgH2, as is evident from figure El Evidently, the 
decomposition is kinetically hindered as it does not occur experimentally. Upon 
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Figure 6. Formation energies of Mg^Ti(-i_^-)H2 as function of composition x for 
ordered fluorite-type (filled red squares) and rutile-type (black stars) structures. The 
lines are added to guide the eye. The open blue squares indicate the formation energies 
of fluorite-type SQSs. 



dehydrogenation the cubic Mg^Ti(i_2,)H2 converts into hep Mg^Ti(i_^). Although the 
latter is metastable with respect to Ti and Mg, its decomposition is also kinetically 
hindered, as thin films can be (de)hydrogenated reversibly. The apparent stability of 
these compounds is quite remarkable. 

The formation energies of the fluorite SQSs are also given in figure [61 Within the 
constraints of the SQS approach, some variation in the relative ordering of Ti and Mg 
atoms is still possible at each composition x. We estimate that this variation gives a 
spread in the energies of ~ 0.1 eV/f.u.. From figure [6] one also observes that the most 
stable SQSs are slightly more stable than the ordered structures, by up to ~ 0.1 eV/f.u. 
for X = 0.875. We have also constructed rutile SQSs of MgJTi(^i_j.) for x = 0.25-0.875. 
However, upon optimization these structures turn out to be unstable and spontaneously 
convert into a fluorite-type structure (i.e. using the conjugate gradient algorithm). 
This would suggest that a fluorite-type structure is stable at least up to a composition 
X = 0.875. As the energies of the two types of structures are quite similar over a range 
of compositions close to the critical composition Xc, the exact value of Xc is difficult to 
determine. Note that here we only consider simple ordered and SQSs models of rutile 
and fiourite-type lattices. Of course, it cannot be excluded that different structural 
models describe the true ground state. However, our structures are judicially chosen so 
that our studies are relevant to understand the experiments; indeed we obtain a good 
agreement with experimental studies. 

The ZPE corrections destabilize the simple hydrides a-MgH2 and a-TiH2, see 
table [21 and one expects this also to be the case in Mg^Ti(i_j,.)H2. The size of the 
ZPE correction in the simple hydrides indicates that it decreases with increasing x 
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from 0.19 eV/f.u. at x = to 0.10 eV/f.u. at x = 1. We suggest that this decrease 
is approximately hnear in x. An exphcit calculation on the SQS of the x = 0.75 
composition gives a ZPE correction of 0.12 eV/f.u., in agreement with this suggestion. 

We now consider to what extend Mg^Ti(i_2.)H2 is suitable as a hydrogen storage 
material. Experimentally it has been observed that the dehydrogenation kinetics 
becomes markedly slower if x > Xq = 0.8 . The results of figure [6] strongly suggest 
that Xo = Xc, i.e. the composition at which the phase transition between fiuorite and 
rutile structures takes place. Lightweight materials require a high content of magnesium, 
but to have a stable fiuorite structure it should not exceed the critical composition 
Xc- We focus upon the composition Mgg 75X10.25112, which should have the fiuorite 
structure. It has gravimetric and volumetric hydrogen densities of 6.2 wt% and 0.135 
kg/L, respectively. Its calculated formation energy is —0.59 eV/f.u. (without ZPE), 
which means that this compound is slightly less stable than a-MgH2 with respect to 
decomposition into the elements. A decreased stability of the hydride should ease the 
dehydrogenation reaction. However, in order that the process is reversible, it is not 
advantageous if the alloy Mg^.Ti(i_2^) dissociates upon dehydrogenation. If we calculate 
the hydrogenation energy according to the reaction 

Mg,.Ti(i_,) + R2i9) ^ Mg,Ti(i_,)H2, (5) 

we find —0.76 eV/f.u. for x = 0.75. This number is in good agreement with the 
experimental value of —0.81 eV/f.u. of [66], obtained if the thin film correction suggested 
there is included. This hydrogenation enthalpy is somewhat lower than that of pure Mg, 
showing that alloying Mg with Ti does not improve the thermodynamics as compared 
to pure Mg. On the other hand, it is also not much worse, and the kinetics of 
(de) hydrogenation of Mgg 75Tio.25H2 is markedly better than that of pure Mg. 

4-3. Electronic structure and optical properties 

The calculated electronic density of states (DOS) of Mg2,Ti(i_^.)H2 in the simple ordered 
structures is shown in figure [7] for fiuorite- and rutile-type structures. Qualitatively 
these DOSs can be interpreted as linear combinations of DOSs of the simple hydrides, 
see figure |2l In the DOSs of the rutile structure one can still distinguish the band 
gap that originates from MgH2. The states just above the gap originate from Ti, as 
discussed in Sec. |3l In particular the DOS of rutile- MgQ g75Tio.i25H2 can be interpreted 
as Ti-doped Q;-MgH2. Enlarging the Ti content then increasingly fills the gap with Ti- 
derived states and for rutile-Mgp 25Tio.75H2, the DOS is dominated by a Ti-derived peak 
at the Fermi level. 

A similar qualitative interpretation of the DOSs also holds for the Mg^Ti(i_^.)Il2 
fiuorite structures. It is somewhat less clear because the band gap in (fluorite-type) 
j3-MgB.2 is smaller than that of (rutile) a-MgH2. Still the DOS of Mgp 875Tio.i25H2 can 
clearly be interpreted as Ti-doped /3-MgIl2, where Ti-derived states appear in the gap. 
Upon increasing the Ti content, the DOS more and more resembles that of /-TiH2 and 
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Figure 7. Electronic density of states (DOS) per formula unit for the fluorite (left 
column) and rutile (right column) Mg^Ti(i_a;)H2 structures. 




Figure 8. Local density of states (LDOS) for fluorite-Mg^Ti(i_j,-)H2, x ~ 0.75 (left 
column) and x — 0.875 (right column). The origin of the energy scale is set at the Fermi 
level. Atomic s, p and d contributions arc shown in red, green and blue, respectively. 
The regions under the s curves are shaded and the d curves are dashed. 
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the gap is completely filled. The latter has implications for the optical properties, as we 
will discuss below. 

In order to be able to identify the character of the electronic states, the 
atom projected DOS of the simple structures is shown in figure [8] for fiuorite-type 
Mg2,Ti(i_2,)H2, X = 0.75, 0.875, which are the most relevant structures and compositions 
for studying the optical properties [15] . It is obvious that the peaks at the Fermi level are 
dominated by states of Ti 3d character, as in /-TiH2, see figure |2l For Mgg 875X10.125112 
these peaks are narrow, reflecting the large distance of 6.62 A between neighbouring Ti 
atoms. One observes the eg-t2g splitting that is typical of d states in a cubic crystal 
field. In this case the latter originates from the cubic coordination of each Ti atom by 
H atoms. Local magnetic moments of 1.5 /iB develop on the Ti atoms and we predict an 
anti-ferromagnetic ordering of these moments at low temperature. The estimated Neel 
temperature is only 40K, however, demonstrating that the moments are atomic like. If 
we increase the Ti content, the states on neighbouring Ti atoms interact more strongly. 
For Mgg 75Tio.25H2 the nearest neighbour distance between Ti atoms is 4.62 A. The Ti 
d derived peaks are much wider, and the local magnetic moments disappear. 

The states in the PDOS at low energy are dominantly H Is derived states, as 
expected. Interestingly, the states around the Fermi level also have a H contribution. 
These contributions are spread out over an energy range that is much wider in 
Mg^Ti(i_^.)H2 than in the simple hydrides MgH2 and TiH2, see figure El The spread 
occurs in H derived bands of s character, as well as in the those of p character. Since 
transitions between such H derived bands give a large contribution to the optical 
response of a hydride [52], one might expect the response to be spread over a large 
energy range. 

We calculate the optical response functions of Mg^Ti(i_^.)H2 as described in Sec. |2j 
The refractive index n{uj) and the extinction coefficient k{uj) are shown in figure [9] for 
the compositions x = 0.5, 0.75, and 0.875. Most relevant for a comparison to the optical 
experiments on thin films are the data for the SQSs. Comparing our data on the SQSs 
and the ordered structures, we observe that one main effect of disorder in the positions 
of the metal atoms, as modelled by a SQS, is smoothing the spectrum. 

In addition, disorder strongly affects the lower end of the spectrum, i.e. u = 0-1.5 
eV, in particular for the Mg-rich compositions x = 0.75, 0.875. The extinction coefficient 
K,{uj) of the SQSs in this frequency range is smaller than that of the ordered structures 
and much closer to its value at high frequency. The large k{u! < 1 eV) in the ordered 
structures is typical for the electron plasma response of metals, which matches the 
DOSs shown in figure [3 As the metallic properties of Mg^Ti(i_a,)H2 originate from the 
Ti atoms only, introducing disorder in the Ti positions can destroy the electron plasma, 
suppressing the large k,{uj < 1 eV) [26]. As can be observed in figure [9l compositions 
with a low Ti content are particularly sensitive to disorder, as the coupling between Ti 
atoms is relatively small. 

Also the structure in the refractive index of the ordered structures at low frequency 
is suppressed in the SQSs and the refractive index at u; > 1 eV becomes close to constant. 
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Figure 9. Index of refraction (top) n{u!) and extinction coefficient k{uj) (bottom) of 
Mg^Ti(i_a;)H2 as function of frequency uj. The black (solid), red (dashed), and blue 
(dotted) lines are for the compositions x — 0.5, 0.75, 0.875, respectively. The left and 
right columns are for the ordered structures and the SQSs, respectively. 

Assuming that one may use bulk optical functions to model the optical properties of 
Mg3,Ti(i_a;)H2, the reflectance (at normal incidence) of a surface is given by the usual 
Fresnel expression R = ((n — l)^ + /t^)/((n + l)^ + /t^). For a metal usually k, ^ n, giving 
rise to a large reflectance. In this case both k and n are rather moderate for frequencies 
in the visible region. Using the data shown in figure [9] one obtains an almost constant 
low refiectance R ~ 0.2 over the whole visible region. Since the extinction coefficient 
is also almost constant in the visible region and it has a quite considerable value, the 
material Mg^Ti(i_a,)H2 has a dark, colourless, i.e. black appearance. 

Using the extinction coefficient of the Mg^Ti(i_2:)H2 SQSs shown in figure M to 
calculate the absorption a{u) = 2k{uj)uj/c, we obtain a good qualitative agreement 
with experimental data [15]. The quantitative difference between the calculated and 
experimental absorption is ~ 30 %, which reflects both the approximations made in the 
calculations, and in the extraction of the experimental values. 

5. Summary and conclusions 

MgH2 has a high hydrogen storage capacity, but it suffers from poor (de)hydrogenation 
kinetics and a high thermodynamic stability, which make it unsuitable as a hydrogen 
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storage material. Adding Ti to Mg gives an alloy hydride Mg^Ti(i_a,)H2 with much faster 
(de)hydrogenation kinetics for x < 0.8. In this paper we have studied the structure and 
stability of Mg^Ti(i_3;)H2, x = 0-1, by first-principles calculations at the level of density 
functional theory. 

As the most stable structures (at ambient conditions) of MgH2 and TiH2 are rutile 
and fluorite, respectively, we model such structures for the alloy hydride Mg^Ti(i_a,)H2. 
A series of simple, ordered fluorite structures is constructed analogous to the high 
pressure MgyTiHie phase, and rutile structures are modelled in MgH2 supercells. All 
cell parameters and atomic positions are optimized and care is taken to converge the 
total energies with respect to computational parameters such as the basis set size and 
the /c-point sampling grid. Upon optimization, all structures stay quite close to either 
fluorite or rutile. To model disorder in the positions of the metal atoms, we construct 
a series of SQSs of Mg^Ti(i„2:)H2 starting from fluorite and rutile structures. Upon 
optimization the fluorite structures are stable, but in compositions with a high Mg 
content, such as MgQ 875X10.125 H2, the structure around the Mg atoms starts to resemble 
that of the high pressure phase /3-MgH2. The SQS rutile structures of Mg3,Ti(i_3.)H2 
are unstable for x < 0.875. Prom the formation energies of the ordered, as well as 
the disordered structures, we give evidence for a fluorite to rutile phase transition at 
a critical composition Xc = 0.8-0.9. This correlates with the experimentally observed 
sharp decrease in (de)hydrogenation rates around this composition. 

Thin fllms of Mg^Ti(i_a,) show a remarkable optical transition from reflecting to 
absorbing upon hydrogcnation, which is reversible upon dehydrogcnation. We calculate 
the optical response of Mg^Ti(i_a,)H2, x = 0.5-0.875, in the ordered and disordered 
structures and show that the black absorbing state is an intrinsic property of this 
material. Mg^Ti(i_a.)H2 has a peak at the Fermi level for all compositions x, which 
is composed of Ti d states. However, the resulting metallic plasma only plays a minor 
role in the optical response in the visible range for compositions with a high Mg content, 
as interband transitions interfere to decrease the extinction coefficient k. Moreover, this 
effect is enhanced by disorder in the positions of the Ti atoms, which easily destroys 
the metallic plasma and suppresses the optical reflection even more. The contributions 
of the H atoms to the bands are spread out over a large energy range as a result of the 
diverse coordination by Ti and Mg atoms in Mg^Ti(i_2.)H2. Interband transitions then 
result in an almost constant optical absorption over a large energy range, causing the 
black appearance of these compounds. 
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